Semi-metal-insulator transition on the surface of a topological insulator with in-plane magnetization 
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A thin film of ferromagnetically ordered material proximate to the surface of a three-dimensional topological 
insulator explicitly breaks the time-reversal symmetry of the surface states. For an out-of-plane ferromagnetic 
order parameter on the surface, parity is also broken, since the Dirac fermions become massive. This leads in 
turn to the generation of a Chern-Simons term by quantum fluctuations. On the other hand, for an in-plane 
magnetization the surface states remain metallic. In this work we study the possibility of spontaneous breaking 
of parity due to a dynamical gap generation on the surface. We find that in the absence of interaction between 
the fermions there is no spontaneous gap generation. In the presence of a local, Hubbard-like, interaction of 
strength g, a gap and a Chern-Simons term are generated for g larger than some critical value, g c , provided the 
number of Dirac fermions, N, is odd. For an even number of Dirac fermions the masses are generated in pairs 
having opposite signs, and no Chern-Simons term is generated. Our result offers a possible explanation to recent 
experiments showing a gap opening even when the topological insulator is proximate to a planar ferromagnet. 
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Due to its unique properties, topological insulators (TI) 
(TJ [2) are likely to play a major role as a component mate- 
rial in different types of hetero structures. For instance, with 
a view towards spintronics applications [3], heterostructures 
involving ferromagnetic (FM) materials or magnetic impuri- 
ties have been studied both theoretically lffl-[T3ll and experi- 
mentally lfT4H20l . Assuming a strong spin-orbit coupling, the 
Lagrangian of the TI surface proximate to a FM thin film in- 
ducing a magnetization n = (n x , n y , cr) on it is given by, 



-£ - ftUyodt - if ' (vpV + z/a) - Jcr]i//, 



(1) 



where if/ = [c^,cj r and iff = ifs^y , with y° = cr z , y l = -icr x , 
and y 2 = icr y . The above Lagrangian has a QED-like form 
in d = 2 + 1 dimension with a vector potential a = (n y , -n x ), 
and no time component for the gauge field. The magnetization 
n(7, r) originates from a proximity effect with a FM thin film 
over the TI surface. Note that time-reversal symmetry is only 
broken if cr ± 0. Thus, an induced planar magnetization does 
not break time-reversal invariance. 

If the electrons on the surface of the topological in are as- 
sumed to interact via a local, Hubbard-like repulsive interac- 
tion, 



Xmt = -|(^rV) 2 , 



(2) 



where g > 0, we can introduce an auxiliary field, ao, via a 
Hubbard- Stratonovich transformation to obtain, 

J 2 

£ = if/[iy°(d t + iJao) - if • (v F V + //a) - J&]iff - — a§. (3) 

^8 

Now the electronic sector is fully gauge-invariant. The only 
term in the Lagrangian that (explicitly) breaks the gauge sym- 
metry is the Gaussian factor in the auxiliary field a^. One 
consequence of this term, quadratic in oq, is that g does not 
renormalize. This is a consequence of the gauge symmetry 
of the fermionic sector and can be easily proved using Ward 
identities |21 ]. An easy way of seeing this without making 



explicit use of the Ward identity is to introduce renormalized 
fields \j/ r = Z" 1/2 ^ and ct r - Z~ l/2 a^ and observe that gauge 
invariance of the fermionic sector implies that the renormal- 
ized exchange coupling is J r = ^Z a J . This implies at the 
same time that terms quadratic in cf which are not gauge in- 
variant should not renormalize, yielding J 2 /g r = Z a J 2 /g and 
consequently g r = g. Therefore, g is a good tuning parame- 
ter in our theory. The fact that g does not renormalize will be 
important in our subsequent analysis. 

Note that we are disregarding the long-range contribution to 
the Coulomb interaction. Such a long-range contribution has 
been shown to be irrelevant in the long wavelength limit in 
theoretical studies of interacting graphene I22H241 and a sim- 
ilar reasoning also applies here. Our Lagrangian corresponds 
to a restricted Thirring model 12511 , in the sense that only the 
zeroth component of the current j^ = ij/y^ij/ appears squared 
in the interaction. 

From the Lagrangian ^ it is readily seen that in the case of 
an out-of-plane magnetization (cr ± 0), the Dirac fermions 
become massive when the symmetry is spontaneously bro- 
ken, in which case the expectation value (cr) ^0. On the 
other hand, in-plane magnetization (cr = 0) implies massless 
Dirac fermions; see Fig. [T] An interesting question to ask is 
whether a mass can also be generated in this case. The an- 
swer to this question is far less obvious than in the case of 
out-of-plane magnetization. This point relates to similar stud- 
ies in (2+l)-dimensional, massless, QED-like theories in di- 
verse condensed matter physics contexts, ranging from high- 
T c cuprates H [26l |27l E3 to interacting graphene ll30H33l . 
However, in the graphene case only the zeroth component of 
the gauge field is present. 

The case of in-plane magnetization is of great experimen- 
tal relevance. Recently a thin film of FM insulator EuS has 
been successfully grown on the surface of Bi 2 Se3 [20], mak- 
ing the surface of Bi2Se3 ferromagnetic by proximity effect. 
The high magnetic anisotropy of the EuS thin film essentially 
induces an in-plane magnetization on the topological surface. 
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(a) Out-of-plane magnetization: gapped Dirac spectrum 




(b) In-plane magnetization: gapless Dirac spectrum 



FIG. 1. (Color online) Schematic comparison between two types 
of magnetization orientation on the surface of a TL (a) In the case 
of out-plane magnetization the electronic spectrum at the surface is 
gapped, (b) For in-plane magnetization the spectrum is gapless. 



In spite of this high magnetic anisotropy, the experiments of 
Wei et al. [20] find a large gap in the Dirac spectrum. In this 
paper we show that such a gap may be spontaneously gener- 
ated in FM/TI heterostructures by a parity symmetry breaking 
mechanism if interaction effects are taken into account. Our 
analysis shows that the long sought evidence for such spon- 
taneous gap generation in QED-like theories may be within 
experimental reach in TIs. 

QED-like theories in two spatial dimensions feature either 
two- or four-component Dirac fermions. Massless QED in 
2+1 dimensions with four-component Dirac fermions has a 
chiral symmetry, which prevents the addition of a mass term 
137) . Indeed, a term mifnfs in the Lagrangian would explicitly 
break the chiral symmetry [37). However, it can be shown 
that a fermion mass is dynamically generated, breaking in this 
way the chiral symmetry spontaneously. This occurs if the 
number of Dirac fermion species, N, is smaller than some crit- 
ical value [36]. In the context of high-r c superconductors this 
mass generation mechanism is associated to the development 
of antiferromagnetic order (26l |28| |29l 134) . In the graphene 
case it is associated to an emergent excitonic gap 1 2 4l [30l432ll . 

The QED-like theory we consider here features two- 
component Dirac fermions and, for this reason, does not have 
a chiral symmetry [ 36 , 37], as 75 -like matrices can only be de- 
fined in representations featuring four-component spinors [21] 
(see, however, the discussion below on the number of Dirac 
fermion species). Hence, the massless case corresponding to 
an in-plane magnetization (<x = 0) has no more internal sym- 
metry than the massive case corresponding to an out-of-plane 
magnetization. On the other hand, a mass term breaks the 
parity symmetry, such that integrating out the fermions gen- 
erates a Chern-Simons (CS) term [38], inducing in this way a 
topological magnetoelectric effect (4). This is what happens 
in the case of an out-of-plane magnetization, where the gen- 



eration of the CS term significantly affects the magnetization 
dynamics [11]. Therefore, it seems that a dynamical mass 
generation in the case of in-plane magnetization would spon- 
taneously break the parity symmetry. For the case of a theory 
containing N two-component massless Dirac fermions, it was 
found in the past using a large N approach that masses are in- 
deed dynamically generated, but half of them acquire a mass 
+m while the other half acquire a mass -m, implying that the 
parity symmetry is actually not broken (39). A recent more 
refined study at finite N [40] confirms this result for N+ = AL, 
where AL denote the number of fermions with mass ±m, re- 
spectively. As a consequence, no CS term is generated, since 
the CS mass is proportional to Y*a=i m al\ m a\ and have no fur- 
ther corrections arising from quantum fluctuations, as required 
by the Coleman-Hill theorem [41]. It is also possible to as- 
sume that some kind of chiral symmetry is broken in the case 
of two-component Dirac fermions by considering instead a 
representation featuring N 12 four-component Dirac fermions 
in a reflection-invariant vacuum state (39l[40). The numerical 
solutions of the Sch winger-Dyson equation (40l in this case 
yield parity violation when (N+ - N-)/(N+ + AL) < 0.4. How- 
ever, this result seems to be in disagreement with a theorem 
by Vafa and Witten [42], which would exclude spontaneous 
parity symmetry breaking in gauge theories of this type. In 
our case the presence of the term proportional to a 2 , in the 
Lagrangian provides a way to spontaneously break parity and 
time-reversal symmetry, since this term breaks both gauge and 
Lorentz symmetries. 

We start our analysis by deriving a purely magnetic action. 
To this end we integrate out N fermionic orbital degrees of 
freedom using an imaginary time formalism and expand the 
logarithm of the fermion determinant up to quadratic order 
in the remaining fields. This leads to the following effective 
action, 



>eff 



1 r cP n r 

2 I On^ U ^{p 2s ^ ~ P»Pv) a »(P) a v(-P) 



-a (p)a (-p) 



(4) 



where p = (co, v F p), the Greek indices run from to 2, and 
U(p) = NJ 2 /(l6\p\) is the vacuum polarization. Here \p\ = 

J co 2 + v 2 p 2 . From Eq. (4 ) we derive the propagator, 



Dfiv(p) = {a^(p)a v (-p)) 



J 2 



p 2 U(p) + 1 



P 2 U(p)\^ y 

PfiPv (/?/Ao + Pv<W 



. (5) 



Interestingly, this result shows that a vector field with a mass 
term only along the temporal direction is not gapped. With an 
isotropic mass term of the form (J 2 1 g)a jl a jl , we would obtain 
instead D^(p) = [p 2 U(p) + J 2 lgY l [^v + (g/J 2 mp)p»Pvl 
which is clearly gapped. As we will see shortly, this difference 
is important, as in our case two of the components of the vec- 
tor field relate to the magnetization and magnetic excitations 
are supposed to be gapless. 



The propagator ^ does not smoothly connect to the 
strongly coupled regime, g — > oo. This is a typical behav- 
ior for massive vector fields |21 ] which is also reflected here, 
although our vector field is only massive along the temporal 
direction. Note, however, that in the strongly coupled regime 
our model reduces to a QED model in 2+1 dimensions with 
two-component Dirac fermions. As mentioned earlier, no CS 
term is generated in this case when the Dirac fermions become 
gapped l39l 

The purely magnetic effective action is finally obtained by 
integrating out ao in the effective action ([4]). This yields, 

S™ = \ J ^ n 0>> [P% ~ F (p)v 2 F PiPj] a,(p)aj(-p) 

= \f ^ n 0>) {[p 2 ~ F(pKv 2 ] »(P) • n(-p) 
+ F(p)v 2 F [p-n(p)][p-n(-p)]}, 



where 



F(p) 



p 2 U(p) + J 2 /g 
v 2 V 2 n(p) + J 2 /g' 



(6) 



(7) 



Going back to real time, the magnetic susceptibility x(u, p) = 
(n+(a),j*)n-(a),j*)), where n+ = n x ± n y , is determined from 
Eq. ([6]) as, 



XfaP) 



16 



NJ 2 Jv 2 p 2 -(a) + id) 2 



2«2 



X 1- 



Ngv A F y 
(co + id) 2 



1 + 



16 



i 



V 2 p 2 - (6o> + /£) 2 



Ng 



v 2 F V 2 



(8) 



From the pole of x(co, p) we infer that the spin- wave velocity 
is identical to the Fermi velocity. Furthermore, by comparing 
with the scaling behavior x((*>, p) ~ [v^p 2 - (co + id) 2 ^ 112 ' 1 for 
a; near v^|p|, yields an anomalous scaling dimension rj = 1. 
This induced anomalous dimension on the topological surface 
is very different from the one of a two-dimensional planar FM 
at T = corresponding to a three-dimensional (d = 2 + 1) XY 
universality class, having rj « 0.04. 

Next, we consider the fermionic propagator. Within an 
imaginary time formalism, the fermion propagator G(p) is 
given in general form by 



/d 3 k 
——^y^G(p - k)D^ v (k)Y v (p, k), 



(9) 



where T v (p, k) is the vertex function. It is understood that 
the Dirac matrices above are the imaginary time counterparts 
of the real time ones defined earlier. They are assumed to 
satisfy the Clifford algebra y^y y + %,% = 26^. In order to 
determine G(p) approximately, we make the decomposition 
G~ l (p) = Z(p)iy fi p fi + H(p) and assume the lowest order form 
for the vertex function, Y^(p,k) = Jy^. Furthermore, we will 
set Z(p) « 1 for G(p) inside the integral in Eq. ([9]). A mass 
will be generated if m = S(0) does not vanish. Note that a 



non-vanishing m implies that (if/ifs) ± 0. Since \jnj/ = n^ - n^ 
mass generation implies also an emergent third component of 
the magnetization. Our strategy will be to make an approxi- 
mation in which H(p) is uniform, H(p) = S(0) = m, and see 
whether there is a solution to Eq. ([9} with m ± 0. We will 
solve Eq. ^ under the assumption that m <*c \p\ «; A, where 
A « NJ 2 is an ultraviolet cutoff. The fermion mass modi- 
fies the vacuum polarization. However, under the assumption 
m <sc \p\ <sc A, the contribution that is even under parity and 
time-reversal, corresponding to the transverse term in Eq. Q, 
remains unchanged. On the other hand, an additional term 
odd under parity is generated. In order to investigate the gap 
equation, we follow Ref . 1 39 ] and assume that N-L fermions 
acquire a positive mass +m, while the remaining L fermions 
acquire a negative mass -m. Thus, the effective action ([4]) for 
the vector field receives the following additional contribution 
odd under parity and time-reversal, 



oodd 
°eff 



AT 



N „ 

2>j 



d 3 p 

(2^ 



U(p)e fivA p A a fi (p)a v (-p). (10) 



This leads in turn to an additional term in the vector field prop- 
agator given by D$\p) = -32^ i (m i /N)e flvA p A /(NJ 2 \p\ 3 ). 
Thus, the following self-consistent equation for m* is obtained, 



1 



16 

~N 



d 3 k 



r d y \ 

J (2^r 



1 



d 3 k 



- 2 — 
m 



) 3 \k\(k 2 + m 2 ) 
N-2L\ C d 3 k 



r d 5 ) 

J (2* 



1*1 



) 3 oj 2 (k 2 + m 2 ) 



N 



r d 3 i 

J (2^r 



1 



) 3 \k\(k 2 + m 2 ) 



d 3 k 



r jPi 



(11) 



t) 3 co 2 (k 2 +rn 2 )' 

The integral in co in the last term of Eq. ( pT) needs to be 
regularized due to the pole at co = 0. To this end, we add a 
small imaginary part co — » co + id and let 6 — > at the end. 

After performing all integrals, the gap equation for i = 
1,. . . , N - L becomes, 



1 



16 



H) 



ln(m/A) + 



4n 



1 + 



1 



while for / = N - L + 1 , . . . , N, we obtain, 



16 



H) 



ln(m/A) + 



gA 
4n 



1 + 



A 



(12) 



(13) 



These gap equations are only compatible with each other if N 
is even, i.e., L = N/2. If we introduce the dimensionless quan- 
tities m = m/A and g = Ag, we obtain for an even number of 
Dirac fermions, 



i*i~ii- 



where 



4;r 2 

7T+1 



81n2 

~n 2 N 



(14) 



(15) 



Thus, we obtain that if N is even a gap is generated provided 
g > g c . In this case there is no generation of CS term. There- 
fore, parity and time-reversal symmetries remain preserved. 
The above result distinguishes itself from the QED case |39l 
due to a complete cancellation of the logarithm. 

The above analysis does not necessarily imply that there 
is no mass generation for an odd number of Dirac fermions. 
Actually, for all N there is in principle the possibility that m* = 
+m for all /. In this case the logarithmic term survives and 
dominates for \m\ «: A over the linear term in \m\. Thus, the 
gap equation is given by Eq. §V2) with L = and where the 
term proportional to \m\/A is neglected. Therefore, we obtain, 



m = exp 



(n+l)N 

64~ 



(g ~ gc) 



(16) 



which is only meaningful for g > g c , otherwise it does not 
decrease with increasing N, which would at large N contradict 
the condition m «: \p\ <$c A, in a situation reminiscent from 
the QED case (39). However, in our case it is possible to 
overcome this difficulty and obtain in addition the generation 



of a CS term. Although Eq. ( 16 ) is in principle valid for any 
N, including the even ones, the solution ( [T4| ) is energetically 
more stable for N even, as it still gives a nonzero gap even in 
the strong-coupling regime g — > oo. 

In conclusion, we have shown that in topological insulators 
with a proximity-induced in-plane magnetization a gap can be 
spontaneously generated by tuning the local electronic inter- 
action above a critical value, leading in this way to a semi- 
metal-insulator transition. We find that when the number of 
Dirac fermions N is even, the masses are generated in pairs 
±\m\ and no CS term is generated, so that parity is overall pre- 
served. On the other hand, for N odd all generated masses are 
equal and positive, so that a CS term is generated. That no CS 
term is generated for N even is physically reasonable, since 
in this case we can change to a representation where there are 
N/2 four-component Dirac spinors, in which case the model 
may be reinterpreted as some model for graphene, a material 
featuring an even number of Dirac cones. Interestingly, in 
such a graphene model the gap generation is associated to a 
mass spectrum containing masses ±\m\, a scenario not consid- 
ered so far in interacting models for graphene where the vector 
field has only the time component [22, 23]. TIs, on the other 
hand, have an odd number of Dirac cones. In this context, 
recent experiments on EuS/Bi2Se3 heterostructures open the 
possibility that the experimentally elusive gap generation in 
QED-like theories may finally be observed in the near future. 
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